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Abstract 



We study the effect of edge contractions on simplicial homology because these contractions have 

turned to be useful in various applications involving topology. It was observed previously that contracting 

O , edges that satisfy the so called link condition preserves homeomorphism in low dimensional complexes, 

.^ and homotopy in general. But, checking the link condition involves computation in all dimensions, and 

hence can be costly, especially in high dimensional complexes. We define a weaker and more local 

^^ condition called the p-link condition for each dimension p, and study its effect on edge contractions. 

,__i We prove the following: (i) For homology groups, edges satisfying the p- and {p — l)-link conditions 

r^ can be contracted without disturbing the p-dimensional homology group, (ii) For relative homology 

r ^ groups, the {p — 1)-, and the {p — 2)-link conditions suffice to guarantee that the contraction does not 

introduce any new class in any of the resulting relative homology groups, though some of the existing 

fj classes can be destroyed. Unfortunately, the surjection in relative homolgy groups does not guarantee 

I— ' that no new relative torsion is created, (iii) For torsions, edges satisfying the p-link condition alone can 

^^ be contracted without creating any new relative torsion and the p-link condition cannot be avoided. The 

►^ results on relative homology and relative torsion are motivated by recent results on computing optimal 

■^ homologous chains, which state that such problems can be solved by linear programming if the complex 

\^ has no relative torsion. Edge contractions that do not introduce new relative torsions, can safely be 

^^ availed in these contexts. 

o 
^' 

O 1 Introduction 

m 

• • The study of edge contractions in the context of graph theory [22], especially in graph minor theory [16] 

. 1^ have resulted into many beautiful results. The extension of edge contractions to simplicial complexes where 

/\^ the structure not only has vertices and edges, but also higher dimensional simplices has also turned out to 

be beneficial for shape representation in graphics and visualization [14, 11] and recently in topological data 
analysis [1]. In this paper, we present several results about edge contractions in simplicial complexes that 
can benefit extraction of topological attributes from a shape or data representation. 

Topological attributes such as the rank of the homology groups, also known as betti numbers, and cycles 
representing the homology classes carry important information about a shape. Naturally, efforts have ensued 
to compute them efficiently in various applications. Examples include computing topological features for 
low dimensional complexes in graphics, visualizations, and sensor networks [3, 7, 20, 23], and for higher 
dimensional ones in data analysis [12]. As the input sizes in these applications grow with the advances in 
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data generation technology, methods to speed up the homology computations become more demanding. For 
example, the Vietoris-Rips complex has been recognized as a versatile data structure for inferring homolog- 
ical attributes from point cloud data [12]. However, because of its inclusive nature, this complex tends to 
contain a large number of simplices and the computation becomes prohibitive when the input point set is 
large. One way to tackle this issue is to use edge contractions [8, 11] that contract edges and collapse other 
simplices as a result. Quite naturally, edge contractions have already been proposed to tame the size of large 
Rips complexes [1]. In a recent work on computing optimal homologous cycles (OHCP), we have shown 
how relative homology and torsions play a role in guaranteeing a polynomial time optimization [9]. To avail 
the benefit of edge contractions in this context we need to understand its effect on relative homology and 
torsion. Motivated by these applications, we make a systematic study of the effects of edge contractions on 
simplicial homology. 

The effects of edge contractions on topology was initially studied by Walkup [21] for 3-manifolds and 
then by Dey et al. [8] for more general domains. They showed that an edge e in a 2-complex or a 3-manifold 
can be contracted while preserving a homeomorphism between the complexes before and after the edge 
contraction if e satisfies a local condition called the link condition. The condition, due to its locality, is easily 
checkable at least for low dimensional complexes. Attali, Lieutier, and Salinas [2] showed that the result can 
be extended to the entire class of finite simplicial complexes if only homotopy instead of homeomorphism 
needs to be preserved. Since homotopy equivalent spaces have isomorphic homology groups, link conditions 
also suffice to preserve homology groups. However, verification of the link condition requires checking it at 
every dimension, which becomes costly for higher dimensional simplicial complexes. 

In this work, we extend the above results in two directions: (i) we study edge contractions not only for 
homology groups, but also for relative homology groups and torsion groups, and (ii) we define a weaker 
and even more local condition called the p-link condition for each dimension p and analyze its effect on 
homology, relative homology, and torsion. Specifically, we prove that an edge contraction cannot destroy 
any homology class in dimension p if the edge satisfies the p-link condition alone. Furthermore, no new class 
is introduced if the edge satisfies the {p — l)-link condition. For relative homology, we show that no new 
homology class in dimension p relative to a particular subcomplex is created if the edge satisfies the (p — 1)- 
and the {p — 2) -link conditions in that subcomplex. This result also implies that no new relative homology 
classes are generated in the contracted complex as long as the contracting edge satisfies the (p — 1)- and 
(p — 2)-link conditions in the original complex. Unfortunately, this surjectivity in relative homology classes 
does not mean surjectivity in relative torsions, that is, the torsion subgroups of relative homology groups. 
Of course, if the edge additionally satisfies the p-link condition in the subcomplex, we have isomorphisms 
in relative homology groups which necessarily mean that no new relative torsion is created. We strengthen 
the condition to require the p-link condition alone if only relative torsion is of interest, which is the case 
for OHCR An example shows that one cannot avoid the p-link condition if one is to guarantee that no new 
relative torsion is introduced. 

Our result on homology preservation under p- and (p — l)-link conditions can be used to compute the 
betti numbers of large complexes after edge contractions. They can also be availed to compute actual rep- 
resentative cycles in a small contracted complex, which can then be pulled back to the original complex 
through a systematic reversal of the contractions. Computations are saved by checking the link conditions 
for only a few dimensions instead of all dimensions. Similarly, our result for relative homology can be used 
to compute the betti numbers of the quotient complexes formed by the original complex relative to a sub- 
complex. Our result on torsion can be availed for computing the shortest cycle in a given homology class. 
This problem, termed the OHCR is known to be NP-hard in general [5]. It has been recently shown that 
the OHCP is solvable in polynomial time when the homology is defined over Z and the simplicial complex 
does not have relative torsion [9]. Similar results hold for related problems such as the optimal bounding 
chain problem (OBCP) [10] and the multiscale simplicial flat norm (MSFN) problem [15]. Therefore, edge 
contractions that preserve the absence of, or even better, eliminate relative torsion should be preferred. Fig- 



ure 3 illustrates a case where an edge contraction does eliminate relative torsion and thus allows computing 
an optimal cycle by linear programming as shown in [9] . We use methods from algebraic topology to prove 
most of our results. Interestingly enough, the result on relative torsion requiring the p-link condition alone 
is proved only with graph theoretic techniques. 

2 Background 

We recall some basic concepts and definitions from algebraic topology relevant to our presentation. Refer 
to standard books, e.g., ones by Munkres [17] for details. 

Definition 2.1. Given a vertex set V, a simplicial complex K = K{V) is a collection of subsets {a C V} 
where cj' C o" is in K if cr € K. A subset u G Kof cardinality p+1 is called a p-i'/m/j/ex. If cr' C cr(cj' C a), 
we call a' a. face {proper face) of a, and a a coface {proper coface) of a' . A map h : K ^ K' between 
two complexes is called simplicial if for every simplex {f i, f2, . . . , Wfc} in i^ {/i(wi), /i(f2), . . . , /i(wfc)} is 
a simplex in K' . 

An oriented simplex a = {vq, wi, • • • , Vp}, also written as v^vi • • -Vp, is an ordered set of vertices. The 
simplices ai with coefficients ai in Z can be added formally creating a chain c = SjajCTj. These chains 
form the chain group Cp. The boundary dpa of a p-simplex o", p > 0, is the {p — l)-chain that adds all 
the {p — 1) -faces of a with orientation taken into consideration. This defines a boundary homomorphism 
dp : Cp ^ Cp_i. The kernel of dp forms the p-cycle group Zp{K) and its image forms the {p — 1) -boundary 
group Bp^i{K). The homology group ^p{K) is the quotient group Zp{K)/Bp{K). Intuitively, a p-cycle 
is a collection of oriented p-simplices whose boundary is zero. It is a non-trivial cycle in Hp, if it is not a 
boundary of a (p + 1) -chain. 

For a finite simplicial complex K, the groups of chains Cp{K), cycles Zp{K), and Hp(i^) are all finitely 
generated abelian groups. By the fundamental theorem of finitely generated abelian groups [17, page 24] 
any such group G can be written as a direct sum of two groups G = F ®T where F = (Z © • • • © Z) and 
T = (Z/ti © • • • © Z/tfc) with tj > 1 and U dividing tj+i. The subgroup T is called the torsion of G. If 
T = 0, we say G is torsion-free. 

A simplicial map h : K ^ K' between two simplicial complexes induces a homomorphism between 
their homology groups which we write as /i* : Hp{K) — ;• Hp{K'). The edge contractions that we deal with 
define such homomorphisms whose properties in relation to various p-link conditions are the focus of our 
work. 

2.1 Link conditions 

We first define edge contractions formally. 

Definition 2.2. Let ab = {a, b} be an edge in a simplicial complex K. An edge contraction of i^ is a 
surjective simplicial map 'jab '■ K ^ K' induced by the vertex map h : V{K) — ;■ V{K') where h is identity 
everywhere except at b for which h{b) = a. 

Authors of [8] investigated when such an edge contraction results in any change in the topology of the 
simplicial complex. They provided a sufficient condition termed the link condition, which guarantees that 
topology is preserved for certain simplicial complexes. This condition has been studied further by Attali, 
Lieutier, and SaUnas [2] recently. 

Definition 2.3. The star of a set X C K, denoted St X, is the set of cofaces of all a ^ X. For a subset 
S of K, the closure of S, denoted S, is the set of simplices in S and all of their faces. Then the link of X, 
denoted LkX, is the set of simplices in St X that do not belong to St X. In the left complex of Figure 1, 
the star of the edge ab consists of ab, abd, abe, abde. Its hnk is d, e, de. 



The following definition of the link condition is taken from the work in [8], and we introduce a weaker 
condition called the p-link condition. 

Definition 2.4. An edge ab ^ K satisfies the link condition in K if and only if Lk a n Lk 6 = Lk ab. It 
satisfies the p-link condition in K, if and only if either (i) p < 0, or (ii) p > and every {p — 1) -simplex ^ 
G Lk a n Lk 6 is also in Lk ab. 

The p-link conditions are weaker than the link condition in the following sense. 

Proposition 2.5. For any edge ab in a simplicial complex K, ab satisfies the link condition if and only if it 
satisfies the p-link condition for all p < dini(i^). 

Proof. We prove both directions by contrapositive. Assume ab does not satisfy the link condition. By the 
definition of star, St a6 = St a n St b. Therefore, St a6 C St a n St b, and Lk a6 C Lk a n Lk b. Therefore, 
there must be a simplex ^ G Lk a n Lk 6 where ^ ^ Lk ab. Let p = dim(^) + 1. Then, ab does not satisfy 
the p-link condition. By the definition of link, there must be a p-simplex ti (^ K where ti = .^ U o, and also 
ap-simplex T2 ^ K where T2 = ^Ub. Therefore, p < dim.{K). 

Now assume ab does not satisfy the p-link condition for some p < dim(K). Then there is some {p — 1)- 
simplex ^ G Lk a n Lk 6 that is not in Lk ab. Therefore, ab does not satisfy the link condition. D 

As an example, consider the two complexes in Figure 1. 
Both complexes have the lower tetrahedron adjoining 
the edge ab. The complex on the left has the top triangle 
missing. The link of the edge ab contains d,e,de in the 
left complex, and it satisfies the 2-link condition but not 
the 1-link condition, ab satisfies both the 1-link and the 
2-link conditions in the complex on the right, since the 
link of ab contains c,d,e, and de. 

3 Homology Preservation 




In this section we study how the p-link conditions for 

edges being contracted affect the homology groups. Not 

surprisingly, the homology classes in dimension p are 

maintained intact if the p- and (p — l)-link conditions 

hold. Specifically, we prove that the p-link condition 

alone implies that no homology class in Hp{K) is destroyed (injectivity) and the (p 

alone implies that no new homology class is created (surjectivity). 




Figure 1: p-link conditions for p = 1, 2. 



l)-link condition 



Tiieorem 3.1. Let ab be an edge in a simplicial complex K and "fab '■ K ^ K' be an edge contraction. 
Then the induced homomorphism jab* cit the homology level has the following properties: 

1- lab* '■ Hp(-ftr) —7- Hp{K') is injective ifab satisfies the p-link condition. 

2- lab* i^ surjective ifab satisfies the (p — l)-link condition. 



To prove the above theorem, we use an intermediate complex K constructed as follows. Let the cone 
f * (T for a vertex v ^ K and a simplex cr G i^ be defined as the closure of the simplex a U {v}. The cone 
V *T ior a subcomplex T C K is the complex v*T = {v*a\a^ T}. We construct K = K U {a*Stb). 
In words, K is constructed out of K by adding simplices that are obtained by coning from a to the closed 



star of b. First, observe that the a * Stbin K can deformation retract to a * Lk 6 taking b to a, and we get 

K'. Therefore, we have the sequence K --^ K ^ K' where i and r are an inclusion and a deformation 
retraction, respectively, and jab = r o i; see Figure 2 for an illustration. At the homology level we have 
the sequence of two homomorphisms where the one on right is an isomorphism induced by a deformation 
retract. 

Hp(K) 'A Hp{k) ^ Hp{K'). 

Since jab* = r^ o i^ at the homology level and r* is an isomorphism, we have that jab* is injective or 
surjective if and only if i^, is. 






Figure 2: Edge contraction as a composition of an inclusion and a retraction. 

Proposition 3.2. i* : \-\p{K) — )• \-\p{K) is injective ifab satisfies the p-link condition, and is surjective ifab 
satisfies the {p — l)-link condition. 

Proof. Since p-dimensional homology is determined only by the skeleton up to dimension p + 1, assume 
that K is only {p + 1) -dimensional. Let S denote the subcomplex a* Sib. The Mayer- Vietoris sequence 



Hp{K n 5) ^H*^ Hp{K) e Hp{S) ^"'J'^ 



Hp{k) ^ Hp.i{K r^ S) 



is exact where i,j,k,i are respective inclusion maps and vr is the connecting homomoiphism [13]. We 
examine the maps j^: ,k^,l^ and prove the required properties for i^, . 

First consider that ab satisfies the p-link condition. We examine the the complex K n S. It is the union 
of St b and the added simplices that are cones of a to the simplices in Lk a n Lk b. None of these added 
simplices can be p-dimensional since otherwise the p-link condition is violated. Any closed star has trivial 
homology, and thus Hp(St b) = 0. Since no simplex of dimension p or more is added to create K n S from 
St b, we still have Hp{K n S) = 0. Therefore, im j* and im A;* are trivial. Because of the exactness of the 
sequence, the map i* — ^* is injective. However, ^* is a zero map since S = a * Stb has Hp(S') = 0. It 
follows that i^ is injective as we are required to prove. 

Next, consider that ab satisfies the (p — l)-link condition. By the same logic as above, there is no 
(p — l)-simplex that can be added to St b to create K n S. Therefore, Hp^i{K n S) = which implies that 
im vr is trivial. Because of the exactness, we have that (z* — £*) is surjective which implies i* is surjective 
as i^: is a zero map. D 



4 Relative Homology Preservation 

In this section we study how the p-link conditions affect the relative homology groups. The motivation 
comes from a recent result on the optimal homologous cycle problem (OHCP) whose efficient solution 
depends on relative torsion and hence on relative homology classes [9] . 



First, we set up some background. Let Lq be a subcomplex of a simplicial complex L. The quotient 
group Cp(L)/Cp(Lo) is called the group of relative p-chains of L modulo Lq and is denoted Cp(L, Lq). The 
boundary operator dp : Cp(L) — ;■ Cp_i (L) and its restriction to Lq induce a homomorphism 

Writing Zp(L, Lq) = ker dp ' ' for relative cycles and Bp_i(L, Lq) = im 9p ' ' for relative boundaries, 
we obtain the relative homology group Hp(L, Lq) = Zp{L, Lo)/Bp(L, Lq). 

A pure simplicial complex of dimension p is formed by a collection of p-simplices and their proper 
faces. We consider relative homology groups Hp(L, Lq) where L C K and Lq C L are pure subcomplexes 
of dimensions (p + 1) and p, respectively. 

To study the effect of edge contraction on relative homology groups, we note that the simplicial map 
7ab : K ^ K' naturally extends to pairs as 'jab ■ {L, Lq) — )■ (L', Lq) where ^ab{L) = L' and ^ab{Lo) = L'q. 
Notice that since L is a subcomplex of K, it may not satisfy some link conditions even if K does. Therefore, 
it is possible apriori that relative homology groups may not be controlled by putting conditions on the edges 
being contracted in K. Nevertheless, we show that for every pair of subcomplexes {L\ Lq) in K' , there is a 
pair of subcomplexes (L, Lq) in K such that ^ab{L, Lq) = (L', Lq) and '^ab* '■ Hp(L, Lq) — )■ Hp(L', Lq) is 
surjective if ah satisfies the {p — 1)-, and the {p — 2)-link conditions in K. This implies that no new relative 
homology classes are created by such an edge contraction. In contrast, however, relative homology classes 
can be killed by edge contractions even if it satisfies all link conditions. We develop these results now. 

The following result known as five lemma in algebraic topology [18, Theorem 5.10] lets us connect our 
result for homology to relative homology groups. 

Theorem 4.1 ([18].). Let h : (L, Lq) — )■ (L', Lq) be a simplicial map. Letfi : Hj(L) — ;■ Hp(L') and gi : 
Hj(Lo) — )■ Hp(Lq) be the homomorphisms induced by h in the absolute homology groups for i = p,p — 1. 
The following statements hold: 

(i) If fp and Qp-i are surjective, and /p-i is injective, then /i* : Hp(L, Lq) — )• Hp(L', Lq) is surjective; 

(ii) If fp and Qp^i are injective, and Qp is surjective, then h^ : Hp(L, Lq) — t- Hp(L', Lq) is injective; 

(Hi) If fi and gi are isomorphisms for i = p,p— 1, then h^ : Hp(L, Lq) — )■ Hp(L', Lq) is an isomorphism. 

If we take the map h to be the restriction of the simplicial map 7^5 to L and its subcomplex Lq, we can 
use Theorem 3.1 to arrive at the following result: 

Theorem 4.2. Let L (^ K be any pure subcomplex of dimension p + 1 and Lq C L be any of its pure 
subcomplexes of dimension p. Let (L', Lq) = ^abi^, Lq). Then, the following hold: 

(i) If ah satisfies the [p — 2)- and (p — l)-link conditions in Lq, 
thenjabt '■ Hp(L,Lo) — )• Hp(L',Lq) is surjective; 

(ii) If ah satisfies the (p — 1)- andp-link conditions in Lq, 
then Jab, '■ Hp(L, Lq) — t- Hp(L', Lq) is injective; 

(Hi) If ah satisfies the (p — 2)- and (p — 1)-, andp-link conditions in Lq, 
then jabt '■ Hp(L, Lq) — t- Hp(L', Lq) is an isomorphism. 



Proof. We prove only (i) from which the proofs for (ii) and (iii) become obvious. If ab satisfies the i-link 
conditions for i = p,p — 2,p — 1 m Lq, then it satisfies the same conditions for L as well. Now apply 
Theorem 3.1 from previous section to conclude the following: 



lab* 
lab* 
lab* 



Hp(L) — )■ Hp(L') is surjective if ab satisfies (p — l)-link condition in Lq 
Hp_i(Lo) —5- Hp_i(Lo) is surjective if ab satisfies {p — 2)-link condition in Lq 
Hp_i(L) — )■ Hp_i(L') is injective if a5 satisfies {p — 1 ) -link condition in Lq 



Now apply Theorem 4.1(i) to finish the proof of (i). D 

It is important to notice that the required link conditions in Theorem 4.2 need to be satisfied in Lq C 
L C K. It is not true that if ab satisfies the i-link condition in K, then it does so in L and Lq. However, for 
i = p — l,p — 2,we have the following observation which allows us to extend the results in Theorem 4.2(i) 
to the case when ab satisfies {p — 1)- and {p — 2)-link conditions in K. 

Proposition 4.3. Let ^ab '■ K — )■ K'. Let [V , Lq) be any pair of subcomplexes of dimensions p + \ and p, 
respectively, where V C K' and Lq C L'. Let L' contain the vertex a. There exists a pair of subcomplexes 
(L, Lq) of dimensions p + 1 andp, respectively, where L (^ K and Lq d L so that jab '■ (L, Lq) — )■ (L', Lq) 
and, for i = p — l,p — 2, ab satisfies the i-link condition for Lq if it does so for K. 

Proof. Consider the preimage of V under jab and take L as its {p + 1) -dimensional skeleton. Similarly take 
Lq as the p-dimensional skeleton of the preimage of Lq. For i = p — l,p — 2, ii there is a simplex a of 
dimension (i — 1), in Lk a n Lk 6 in Lq, then the simplex ab * a has to be present in Lq since this simplex 
is in K because ab satisfies the i-link condition in K. Therefore, ab satisfies the same link condition in Lq 
as well. n 

Notice that the above observation does not include the p-link condition. Since the p-link condition may 
require a (p + 1) -simplex, the edge ab may not satisfy the p-link condition in the p-dimensional complex 
Lq, even if it does so in K. This is why we cannot extend Theorem 4.2(ii) and (iii) in the result below 
which relates the surjectivity of the relative homology groups with the link conditions in the original input 
complex K instead of a subcomplex. The proof of this theorem follows from applying Proposition 4.3 to 
Theorem 4.2(i). 

Theorem 4.4. Let L' C K' be any pure subcomplex of dimension p+1 and Lq C L' be any of its pure 
subcomplex of dimension p, where K' = jab{K). There exists a pair of pure subcomplexes L, Lq of dimen- 
sions p +1 and p, respectively with L C K and Lq C L so that: 

Ifab satisfies (p — 2)- and {p — l)-link conditions in K, then jab, ■ Hp(L, Lq) — t- Hp(L', Lq) is surjective. 

4.1 Implications of Theorem 4.4 and Theorem 4.2 

Surjectivity. Notice that Theorem 4.4 is sufficient to claim that edge contractions cannot create any new 
class in relative homology groups if the edge ab satisfies only (p — 1)- and (p — 2)-link conditions in K. 
However, it may happen that the sujectivity does not respect torsions, that is, the preimage of a torsion 
subgroup in Hp(L', Lq) may not have torsion in Hp(L, Lq). This would mean that (p — 1)- and (p — 2)- 
link conditions are not enough to guarantee that torsion subgroups have sujection only from the torsion 
subgroups of the source space. We are interested in link conditions that guarantee that no new torsion class 
is created by an edge contraction because of its connection to the problem of OHCP as discussed in the next 
section. 



We present an example where a new relative torsion indeed appears in the contracted complex even 
though the edge satisfies the {p — 1)- and {p — 2)-link conditions. Consider a sequence of triangles forming 
a Mobius strip. It is known that Mobius strip has a torsion in Hi relative to its boundary. Now remove 
one triangle, say abc, and call the new complex M. Assume that the edge ab was on the boundary of M. 
The complex M, which is a Mobius strip with one triangle removed, does not have any relative torsion. 
However, if we contract the edge ab, we eliminate the hole created by the removal of abc. The resulting 
complex now is a Mobius strip and hence has a relative torsion. Notice that, for p = 1, ab satisfied {p — 1)- 
and {p — 2)-link conditions vacuously though it did not satisfy the p-link condition. 

The isomorphism result in Theorem 4.2(iii) guarantees that no new relative torsion appears if ab satisfies 
all three link conditions, namely p-, {p — 1)-, and {p — 2)-link conditions in Lq which cannot be guaranteed 
by requiring ab to satisfy them in K. 

In section 6, by a graph theoretic approach, we show that the p-link condition in K alone is sufficient to 
prevent the appearance of new relative torsions. The above example illustrates that the p-link condition is 
necessary for guaranteeing surjectivity in relative torsions. 

Injectivity. Unlike surjection, the injectivity implied by Theorem 4.2(ii) cannot be used to claim that no 
relative homology class will be killed by an edge ab satisfying the p- and {p — l)-link conditions in K. 
The reason is that for injectivity of 7^5, we choose the pair (L, Lq) in K first and then consider its image 
under jab- So, even if ab satisfies p- and {p — l)-link conditions in K, it may not do so in Lq and hence 
Theorem 4.2(ii) may not be applied. In fact, contracting edges satisfying all link conditions in K can indeed 
kill a relative torsion. Figure 3 illustrates such an example. The resulting complex after contracting edge 
ab in the 2-complex on left is shown on right. The left complex minus the triangles abd, abm and the edge 
ab forms a 15-triangle Mobius strip with a self-intersection at vertex d. This strip relative to its boundary 
results in a relative torsion. On the right, the self intersection expands to the edge ad, which causes the 
relative torsion to disappear. 

This example illustrates that an edge contraction may destroy a relative torsion though no new relative 
torsion is generated thanks to Theorem 4.4(iii). This is of course good news for OHCP since such edge 
contractions can only take away obstructions to their solutions and not introduce new ones. We elaborate 
on this statement now, and describe a concrete example illustrating the benefits of such edge contractions to 
the efficient solution of OHCP instances. 

4.2 Relative homology preservation and the OHCP 

Given an oriented simplicial complex K of dimension d, and a natural number p, I < p < d, the p- 
boundary matrix of K, denoted [dp] , is a matrix containing exactly one column j for each p-simplex a in 
K, and exactly one row i for each (p — 1) -simplex r in K. If r is not a face of a, then the entry in row i 
and column j is 0. If r is a face of a, which we denote by r ^ o", then this entry is 1 if the orientation of r 
agrees with the orientation induced by a on r, and — 1 otherwise. 

Note that given a simplicial complex K and a choice of orientations for all simplices in K, [dp] is unique 
under row and column permutations, and will generally be referenced as a single matrix. We also need the 
notion of total unimodularity. A matrix A is totally unimodular, or TU, if the determinant of each square 
submatrix of A is either 0, 1, or —1. An immediately necessary condition for A to be TU is that each 
aij G {0, ±1}. The importance of TU matrices for integer programming is well known - see, for instance, 
the book by Schrijver [19, Chapters 19-21]. In particular, it is known that the integer linear program 

mill f X subject to ylx = b, x > and x € Z" 

for A e Z"*^", b G Z" can always, i.e., for every f G M", be solved in polynomial time by solving its linear 




Figure 3: Contracting ab destroys a relative homology class even though it satisfies the link condition. In 
the left complex, a Mobius strip is formed by the triangles with orientations shown in black. This Mobius 
strip does not exist in the right complex, since it self-intersects along edge ad here. The input chain for the 
OHCP instance is edge ef (shown in purple). Candidates for optimal homologous chains are shown in green 
and cyan on the left complex, and in green on the right complex. See Section 4.3 for more details. 

programming relaxation (obtained by ignoring x G Z") if and only if A is totally unimodular. The main 
motivation for our discussion of totally unimodular matrices is the following result in [9]. 

Theorem 4.5. For a finite simplicial complex K of dimension greater than p, the boundary matrix [i9p+i] 
is totally unimodular if and only if Hp (L, Lq) is torsion-free, for all pure subcomplexes Lq, L in K of 
dimensions p and p+ 1 respectively, where Lq C L. 

Following this result, instances of OHCP and related problems (with homology over Z) can be solved 
in polynomial time when the simplicial complex K is relative torsion-free. 



4.3 An example where edge contraction helps to solve the OHCP 

We present a small example which illustrates the effectiveness of edge contractions on efficient solutions of 
OHCP. Consider the simplicial complexes K on the left and K' on the right in Figure 3. We obtain K' from 
K by contracting the edge ab, which satisfies the link condition. 

We consider the following OHCP instance on K, and equivalently on K'. All edges in red and the edge 
ef in purple have weights of 1 each. All edges in green and the thinner edges in black have weights of 0.05 
each. The two edges bd and bm in K, drawn in cyan, have weights of 0. 10 each. Orientations of all triangles 
and the edge ef are shown. Remaining edges could be oriented arbitrarily. The input 1-chain c consists of 
the single edge ef with multiplier 1. 



The unique optimal solution to the OHCP LP is the chain x consisting of the 15 edges shown in green 
and cyan in K, with multipliers of ±0.5, resulting in a total weight of 0.425. x is homologous to c over 
M, as their difference is the boundary of the 2-chain with multipliers of —0.5 for each of the 15 triangles in 
the Mobius strip, whose orientations are shown using black arrows (triangles abd and ahm have multipliers 
of 0). But the unique optimal homologous chain sought here is the 1-chain x' shown in green in the right 
complex (this chain is identical in K and K'). Consisting of the 9 edges ak, km, ma, ec, en, nl, Ij, jd, df 
with multipliers ±1 each, x' has a total weight of 0.45. But to obtain this solution, we have to solve the 
OHCP model as an integer program, instead of as a linear program. 

Notice that K' obtained from K by contracting edge ab is free of relative torsion, and the weights of 
each edge in K' is the same as its weight in K. The unique optimal solution to the OHCP LP on K' for the 
identical input chain is x' itself, x' is homologous to c as defined by the 2-chain consisting of the triangles 
whose orientations are shown using blue arrows, each with a multiplier of —1. Hence we are able to solve 
the original OHCP instance using linear rather than integer programming, after contracting one edge. 

We show that satisfying the p-link condition alone is sufficient to guarantee that no new relative torsion 
is introduced. Instead of proving the result directly, we show that the p-link condition for the edge being 
contracted preserves the total unimodularity of the boundary matrix, which in turn guarantees the absence of 
relative torsion thanks to Theorem 4.5. In contrast to our earlier approach, we use results from graph theory 
to arrive at this result. 

5 Bipartite Graphs and Boundary Matrices 

Given a matrix A with entries in {0, ±1}, we associate with it a weighted, undirected, bipartite graph 
G = {ViUV2,E) where each row of A corresponds to a node in Vi and each column of A corresponds to a 
node in V2 [4, 6]. Each nonzero entry aij in A is associated with an edge connecting the nodes in Vi and V2 
corresponding to row i and column j, respectively, with a weight equal to a-ij. We call G the bipartite graph 
representation of A. 

Some definitions from graph theory are central to our discussion. For a subgraph S* of G containing a 
vertex v, we denote the number of edges in S incident to v as the degree of v in S, or deg5(f ). A cycle 
y is a connected subgraph of G where for each vertex v ^ Y, degy(t') = 2. We call a subgraph C of G 
a circuit if for each vertex w G C, degc(f ) = mod 2 . By this definition, it is possible that v ^ G, but 
degc(t;) = 0. However, we consider two circuits as equivalent if they contain the same edge sets with the 
same weights for each edge. This means, for example, that a circuit C containing only vertices of degree 
in C is equivalent to the empty subgraph. 

Definition 5.1. A circuit G in the weighted graph G is b-even if the sum of the weights of the edges in G 
is mod 4, and b-odd if the sum of the weights of the edges is 2 mod 4. The quality of G being b-even, 
b-odd, or neither is called the b-parity of G. 

This definition is equivalent to the definition of even and odd cycles given by Conforti and Rao [6]. 
Given a simplicial complex K, consider the p-graph Gp{K) of K constructed as follows. Each p- and 
(p — l)-simplex a £ K provides a dual vertex a* in Gp{K). We call the vertex a* a p- or {p — l)-vertex if 
aisap- or (p — l)-simplex, respectively. There is an edge (t*t* in Gp{K) if and only if the {p — l)-simplex 
a £ K has the p-simplex r as a coface. The weight of a*T* is 1 or — 1 depending on whether a and r match 
in orientation or not, respectively. It is evident that Gp{K) is a weighted bipartite graph whose adjacency 
matrix is given by [dp] . This means the following proposition is almost immediate. 

Proposition 5.2. Reversing the orientation of any collection ofp- and (p — l)-simplices in K does not alter 
the b-parity of any circuit in Gp{K). 

10 



Let (7 be a circuit in a graph G. A chord of C is a single edge of G not in C whose both end points are 
vertices in C. If C has no chord, it is called chordless, and to say G is induced is equivalent to saying G is 
chordless. Using the terminology given above, we now state without proof an important result presented by 
Conforti and Rao [6], who were in turn motivated by the results of Camion [4]. 

Theorem 5.3. For a matrix A with entries in {0, iti} and its bipartite graph representation G, A is totally 
unimodular if and only ifG contains no chordless b-odd circuit. 

By combining Theorem 4.5 and Theorem 5.3, the following corollary is immediate. 

Corollary 5.4. For a finite simplicial complex K of dimension greater than p, the following results are 
equivalent. 

1. Hp (L, Lq) is torsion-free for all pure subcomplexes Lq, Lin K of dimensions p and p+1, respectively, 
where Lq C L. 

2. The boundary matrix [Sp+i] is totally unimodular 

3. The bipartite graph representation Gp+i{K) of [dp+i] contains no chordless b-odd circuit. 

5.1 Simplices and edge contraction 

Recall that each simplex in K is mapped by an edge contraction 'jab : K ^ K' to a. simplex in K'. We 
categorize simplices into three cases based on how they get mapped by 7^5. These cases are defined relative 
to the specific edge ab being contracted. We illustrate these cases in Figure 4, and introduce several related 
definitions below. 

1. For each pair of simplices a, a' G K where a ^ a,b ^ a', and a = {a' \ {b}) U {a}, we have 
7ab{cr) = jab{cr') = a. Then a and a' are mirror simplices, and we say they are mirrors of each other. 
Their duals are mirror vertices, and we say these vertices are mirrors of each other. 

In Figure 4, a is the mirror of b, ad is the mirror of bd, ae is the mirror of be, and ade is the mirror of 
bde. Similarly, dual vertex ui is the mirror of U2, vi is the mirror of V2, and wi is the mirror of W2. 

2. a ^ K is collapsing if a, 6 G a. Its dual a* is a collapsing vertex. Note if a is p-dimensional, then a 
has exactly one pair of (p — l)-faces that are mirrors of each other. Then 7a6(o') = r, where r is the 
unique mirror (p — l)-face of a containing a. 

In Figure 4, ab, abd, abe, and abde are all collapsing simplices. V3, w^, and s are collapsing vertices. 

3. a €^ K is injective if neither of the above cases applies. We have 7afe(o') = cr', and 7~fe^(o"') = a. If 
b ^ a, then a' = a.lib £ a, then a' = {a\ {b}) U {a}. 

We state a few more definitions related to the Mirror and Collapsing cases. 

Definition 5.5. Let ai € Kbea.{p-\- l)-simplex with mirror cj2 G K. Note that cri and 0-2 have exactly one 
common p-face ^ that is an injective simplex. Each other p-face ri ^ cJi is a mirror of some r2 ^ o"2, and 
vice-versa. In any subgraph S of the dual graph Gp+i of K that contains both edges (t^Ti and 0-2 r|, these 
two edges are mirror edges, and are mirrors of each other. In Figure 4, assuming S is the entire 2-graph, tui 
is the mirror of tu2, uivi is the mirror of U2V2, and uiwi is the mirror of U2W2. 

Definition 5.6. For any p-simplices ri and T2 that are mirrors of each other in K, the two unique edges of 
the dual graph Gp+i of K that directly connect ri and T2 each to their common collapsing (p + l)-coface 
a are together called the mirror connection between n and T2, and also between r* and T2. There are two 
mirror connections in Figure 4. The mirror connection between vi and V2 are the two edges v^vi and v^V2, 
and the mirror connection between wi and W2 are W3W1 and W3W2. 
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Figure 4: The 2-graph of a complex, and the result after edge contraction 7a(,. 

Definition 5.7. Any edges incident to two vertices in the dual graph Gp+i that are both collapsing are called 
collapsing edges. In Figure 4, sv^ and sw^ are collapsing edges. 

6 Relative Torsion- Aware Edge Contractions 

We now present the theorem which states that contracting an edge in a simplicial complex satisfying the 
p-link condition does not create any new relative torsion in dimension p. In other words, if the {p + 1)- 
boundary matrix of a simplicial complex is totally unimodular to start with, such an edge contraction will 
preserve this property. Theorem 6.1 below states the same thanks to Corollary 5.4. 

Theorem 6.1. Let K' = ^ab{K) where ah satisfies the p-link condition. Let ^ be the set of circuits in graph 
G = Gpj^i{K), and ^' be the set of circuits in G' = Gp^i{K'). Let ^x, <^ ^ be the set of all circuits 
of G that contain a collapsing p-vertex. Let ^m ^ ^ be the set of all circuits of G that contain a pair of 
(p + 1) -vertices that are mirrors of each other There exists a function / : ^ \ {^m U '^l) — ^ ^' with the 
following properties. 

PI. / is surjective. 

P2. / preserves b-parity. 

P3. For each G in the domain of f, ifG has a chord then so does f{C). 

Proof. Before providing details of the proof, we refer the reader to Figure 4 to visualize '^a/ and 'rfi for 
p = 1. ''i^L is the set of all circuits that contain edges incident to the vertex s, which is the dual of the only 
collapsing p-simplex in this complex. One such circuit is {sv3, V3V1, viui, uiwi, wiw^, wss}. ^m is the 
set of circuits which contain an edge incident to ui, and also an edge incident to U2. One such circuit is 
{uit, tu2, U2V2, V2V3, V3V1, viui}. Though there are several circuits in this graph, none are in the domain of 
/ as all of them belong to either ^m or "^x. 

We break the proof up into several parts. We first define the function / in Section 6.1, and then show 
that / is surjective (Property PI) in Section 6.1.1. Finally, we prove that / preserves 6-parity (Property P2) 
in Section 6.1.2 and that / preserves chords of circuits (Property P3) in Section 6.1.3. 
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6.1 Definition of / 

We define the function / in the following way. For any C € Dom(/), let L = {a\a* G C}. For each 
vertex a* of C, f{a*) = jabio'T- The function / maps each edge a*T* £ C to the edge Jabio-)* 1ab(,T)* 
if the edge is not a mirror connection, otherwise simply to the vertex 7afe(<7)* = 7a6(''")*- Since 7^;, is well 
defined, for any C G Dom(/), if /(C) = X and /(C) = Y, then X = Y. 

To show /(C) G '^', first note that if all simplices of L are injective, this is almost trivially so since in 
these cases, 'jab is injective, preserves the dimension of simplices, and r G L is a face of cr G L if and only 
if labir) is a face of 7afe(o-). 

If L does not contain any collapsing simplices, but contains a pair of mirror simplices ri and T2, then 
since C ^ 't^m, ti and T2 are p-simplices. Furthermore, the proper cofaces of either ri or T2 are injective 
because the only coface which cannot be such is a collapsing simplex. Let r = 'jabi'Ti) = 7afe(T2)- Since 
7ab maps the proper cofaces of ri and T2 injectively, degft(^\{T*) = degc(rj*) + degc(r|). Since both 
operands of this sum are even, the sum is even. Therefore, /(C) G ^'. 

Now assume that L contains a collapsing simplex a. Since C ^ "^l, a is a (p+ 1) -simplex. Furthermore, 
since all but two p-faces ri and T2 of a are also collapsing, a* must be of degree 2 in C. Therefore, C must 
contain the mirror connections T*a* and t|cj*. Note that a is the only {p + 1) -simplex that is a coface of 
both Ti and T2. The edge contraction 7^5 maps other proper cofaces of ri and T2 in L injectively. Both 
mirror connections r^a* and T2cr* are mapped by / to the vertex r* where r = 7afe(Ti) = 7a6(''"2) = 
7a6(o-). Therefore, degj(c)(r*) = deg(;(rf ) + degc(r2*) - 2. Again, since degcirl) + deg(^(T2*) is even, 
deg/(c)(r*) is even. Hence, /(C) G ^'. 

6.1.1 Surjectivity of / 

We have shown that / maps each element in its domain to an element in its codomain. To show that / is 
surjective (Property PI), we define a function that extends the domain of / to all subgraphs of G. Let S^ be 
the set of all subgraphs of G, and J^' be the set of all subgraphs of G'. Define the function g : J^ ^ y 
as the extension of / to this domain. Hence it is possible that some vertex r* G Dom((7) is the dual of a 
collapsing p-simplex r. In this case we specify that g maps r* and all edges of C incident to r* , which all 
must be collapsing edges, to the empty subgraph. Since 7^5 is surjective, so is g. To show / is surjective, we 
will take an arbitrary C G '^' and show that its preimage under g, which we will denote 5^c^ must contain 
an element of Dom(/). 

If 5^c contains an element S ^ Dom(/), there are three type of edges that we may remove or add to S 
without affecting g(S), as detailed below. 

— We may remove or add any collapsing edge. 

— We may remove or add any edge that is a mirror connection in G. 

— We may remove either one, but not both, of any two edges that are mirrors of each other. For any edge 
in S that is a mirror in G but not in S, we may add its mirror to S. 

Using these modifications, we describe a procedure to construct a C/ G Dom(/) n S^c (see Figure 5). 
Figures 6, 7, and 8 illustrate the steps of this procedure. Graph vertices are labeled with the vertex labels of 
the dual simplices. In each figure, graph edges and vertices eliminated by the step shown are highlighted in 
red, and edges and vertices added by the step are highlighted in green. 

We provide some details of the steps in the procedure now. Denote the subgraph after Step II as 82- Since 
we have removed any collapsing edges, the only edges remaining in 5*2 that g does not map injectively are 
mirror connections in G. 

Since G' is a circuit, the only vertices that are of odd degree in S2 are vertices incident to a mirror 
connection in G. 
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Procedure Construct Cj 

Input : S e yc\ Dom(/) in graph G such that g{S) G '^'. 
Output : C/ G ^c n Dom(/) such that g{Cf) = g{S). 

I. Remove any collapsing edges. 

II. For any edge a*T* G S that is a mirror in G, where a adjoins the contracting vertex b, replace it 
with its mirror in G if this mirror is currently absent in S, or remove a*T* otherwise. 

III. For each pair of vertices of odd degree that are mirrors of each other, negate the mirror connections 
connecting them, i.e., remove them if they are in the current subgraph, or add them if they are not. 

IV. For each pair of vertices of odd degree incident to a common single edge in G, negate this edge. 



Figure 5: Procedure to construct a circuit in Dom(/) starting with a subgraph not in Dom(/). 

For any pair of mirror connections connecting two mirror vertices in G, there are three vertices involved 
- two p- vertices rj* and t| that are mirrors of each other in G, and a collapsing (p + 1) -vertex a*. degg^{a*) 
is the number of edges of the mirror connection in 5*2 , and these are the edges g removes by mapping them to 
the vertex f = ^(tj*) = g(T|) = g{cr*). Because C is a circuit, and edges that are not a mirror connection 
are mapped injectively by g, degg^{a*) is odd if and only if degg^{T*) + degg^{T2) is odd. Therefore, 
among the three vertices Ti,T2,t*, an even number of them must be of odd degree in 82- The graph on the 
right in Figure 7 illustrates this situation. Note that the vertices ahi* and bhi* have degrees 3 and 1 in S2, 
respectively, and the vertex abhi* has degree 2 in 5*2. 

If T*,T2, cr* are all of even degree, no further action is needed. If this is not the case, but a* is of even 
degree, then Step III will make r* and t| of even degree, and a* will still be of even degree. Notice that 
since ab satisfies the p-link condition, the mirror connections T^a* and t|(T* must exist in G. If a* is of odd 





Figure 6: A subgraph S of the 3-graph where g{S) = C" G '^' (left) and changes in Step I (right). 
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degree, then Step IV will convert a* and whichever vertex between rj* and r| is of odd degree to vertices of 
even degree. 

Steps III and IV ensure that every vertex in the subgraph have even degrees. The resulting subgraph Cj 
is a circuit in J^c that is not in 't^m because we have removed at least one of the edges adjoining vertices 
that are mirror to each other in Step II, and not in 'tfi because we have removed all collapsing edges in Step 
I. Therefore, Cf e Dom(/) n ^c- 

6.1.2 Preservation of 6-Parity 

To show / preserves 5-parity (Property P2), let C G Dom(/), C = f{C), and L = {a\a* G C}. We 
need to analyze only those vertices a* of C where a is collapsing, or is a mirror. Since C G Dom(/), we 
may restrict our attention to mirror p- vertices in C, or collapsing {p + l)-vertices in C. In fact, we may 
further restrict our attention to edges incident to p-vertices rj* and T2 that are mirrors and in C, because 
C G Dom(/) implies any edges incident to a collapsing {p + l)-vertex is also incident to a p-vertex that is 
a mirror in C. 

We will divide the possibilities into four cases based on two criteria - either the mirror connections 
connecting r* and T2 are both in C or both not in C, and either the orientations of ri and T2 are consistent, or 
they are not. Recall that if ri and T2 have a unique common {p— l)-face ^, ri and T2 are consistently oriented 
if they induce opposite orientations on ^. We assume that unless otherwise specified, the orientations of 
simplices in K' are chosen so that for each edge e in C, the weights of e and g{e) are the same. 

If neither edge of the mirror connection is in C, / maps all edges of C injectively to edges of C'. If 
the orientations of ri and T2 are not consistent, then we may choose the ordering of vertices of r = 7afe(Ti) 
that determines its orientation to be the same as the ordering of the vertices of ri. Then each edge e of C" 
will have the same weight as g~^{e) in C. If the orientations of ti and T2 are consistent, then we may still 
choose the order of vertices of r in the same way as described above, but now each edge e incident to r| 
will have the opposite weight as g{e). The sum of the weights of edges incident to r| and the sum of the 





Figure 7: Changes in Step II (left) and in Step III (right). 
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Figure 8: Changes in Step IV, and Cj. 

weights of their corresponding edges in C must differ by 2 for each such edge e. Since C is a circuit, there 
must be an even number of such edges. Therefore, this difference is mod 4, and hence the 5-parity of C 
is the same as that of C. We now discuss the remaining case when both mirror connections are in C and the 
orientations of ri and T2 are not consistent. We illustrate both cases for p = 1 in Figure 9. 

Consider the case when both mirror connections are in C and the orientations of ri and T2 are not 
consistent. Notice that the orientation of their common coface a must agree with the orientation of one of 
these mirrors and disagree with the orientation of the other mirror. Hence, the sum of the weights of the 
mirror connections must be 0. The function / maps all other edges of C injectively to edges of C", and we 
may choose the orientation of r as before. Therefore, the sum of the weights of edges in C" equals the sum 
of the weights of edges in C. 

If the orientations of ti and T2 are consistent, then the sum of the weights of the mirror connections must 
be 2 mod 4. However, the number of edges mapped injectively by g that are incident to r| must be odd. 
Hence, the difference between the sum of weights of these edges and their corresponding edges in C" must 
also be 2 mod 4. It follows that the overall difference between the weights of edges in C and C is mod 4, 
and the 6-parity of C' equals that of C. Thus, in all cases, / preserves 6-parity. Figure 9 illustrates simple 
examples of the cases when the mirror connection is in C for a 2-graph. 

6.1.3 Preservation of Chords 

To show property P3, notice that each edge of G is contained in a (p + l)-simplex in K. If C € Dom(/) 
has a chord h, and h adjoins a (p + l)-vertex a*, then there must be at least three p-vertices in C - at 
least two directly connected to a* by edges in C, and also the p-vertex adjoining h. Since each collapsing 
{p+ l)-simplex has exactly two p-faces that are not collapsing, a cannot be a collapsing simplex. Otherwise, 

If L represents the set of all simplices in K whose duals are in C, a cannot have its mirror in L as 
C 'ta'M- Therefore, a is injective in L. Figure 10 illustrates a chord in a 2-graph. Since jab niaps a 
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Mirror connection in circuit, orientations not consistent 




sum = 



Mirror connection in circuit, orientations consistent 




sum = 2 diff. = 2 

Figure 9: Examples of weights of edges of the 2-graph before (left) and after (right) edge contraction jab- 



injectively, any edge of G in o" is in C if and only if its image under g is in /(C). Therefore, g{h) is not 
in /(C), and for the vertex a* of h, g{a*) is in /(C). To show for the other vertex r* of h that g{T*) is in 
f{C), this could only be not true if all edges of C incident to r* were not mapped to edges by g. Since C 
is a circuit, r* is incident to at least two edges. But the only p- vertices incident to more than one edge in G 
not mapped to an edge by g are dual to collapsing p-simplices, which cannot hold for r* here since C ^ '^l- 

D 




Figure 10: A chord h of cycle C, and C = /(C). 



7 Discussion 

We have presented several results that connect the local p-link conditions to preserving various homolog- 
ical classes and torsions under edge contractions. We used both graph theoretic and algebraic topological 
techniques to arrive at these results. Our results on homology and relative homology groups may be used 
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to accelerate the computations of these topological structures efficiently using edge contractions. More im- 
portantly, we have laid down almost a complete picture of the relationship of edge contractions in regards to 
simplicial homology. It is not hard to see that the conditions in Theorem 3.1 are not only sufficient, but are 
necessary for absolute homology groups preservation. For relative homology, the scenario becomes more 
subtle. The injectivity and hence isomorphism cannot be guaranteed with any link condition in the original 
complex as our example in Figure 3 shows. The link conditions in Theorem 4.4 are sufficient for surjec- 
tivity; but are they necessary? This remains open. For relative torsions. Theorem 6.1 and our example of 
punctured Mobius strip in Section 4. 1 show that the p-link condition is sufficient and sometimes necessary 
for preventing new torsions, and our example in Figure 3 shows that no link condition can ensure preserving 
existing relative torsions. 

An open question is whether there exist local conditions that also preserve the total unimodularity prop- 
erty of the complex. We know that total unimodularity cannot be destroyed by an edge contraction if the 
appropriate link condition holds. But, we do not know how to preserve the absence of total unimodularity 
with local conditions. In the Section 4.3, we discuss a concrete example where linear programming fails to 
provide an optimal homologous chain, but an edge contraction makes it amenable to such computation by 
eliminating relative torsion. An immediately relevant question is how to approximate the optimal chain in 
the original complex using the one in a reduced complex. 
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